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ORBITAL CHANGES OF THE GASEOUS RING ARDUJ’D STARS 

Su-Shu HuajDg 

Department of Astronomy, Northwestern University 


ABSTRACT 


In order to understand the seemingly erratic V/R variations and radial 
velocity curves of emission edges and central absorption of some Be stars, 
we have advanced, for the changes of structure of the gaseous ring around 
these stars, a theory that is based on the interaction between the existing 
ring and the newly ejected matter from the star. It shows that the struc- 
tural change of the ring is completely controlled by the angular-momentum 
input factor and the dissipation factor. In light of this understanding, 
we have gone on to interpret the observed results of A Mon and 7T Aqr. 
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I. INTRODUCTION 


Be stars differ from B stars mainljk by the presence of emission^ lines 
which are generally agreed to come from the rotating gaseous ring around the 
central star (Struve 1931). Therefore to understand Be stars it is Important 
to know the formation and evolution of these rings. In previous papers (Huang 
1937b, 1975) a scenario has been proposed to show how exactly the ring is 
formed as a combined result of rapid rotation of the central star and the 
radiation pressure that exerts on the matter in the equatorial region. The 
rapid rotation provides the angular momentum while the radiation pressure 
serves the purpose of lifting up the matter from the stellar surface. Unaware 
of Huang's papers, Massa (1975) has independently advanced some of Huang's 
arguments. 

That the radiation pressure plays a critical role in the formation of 
the gaseous envelope around Be stars was pointed out by Gerasimovic (1934) 
a long time ago, but Huang's idea of the existence of multi-substrata in the 
gaseous medium around the star had gone a step further than Gerasimovic 's. 
Actually it is easy to see the necessity of introducing multi-substrata if 
one considers the fact that the specific angular momentum (i.e. angular momentum 
per unit mass) of a particle in Keplerien motion around the star is proportional 
to the square root of semi -major axis of the orbit. So even if the star ro- 
tates at the break-up velocity at its equator, the ejected matter cannot all 
be converted into a rotating ring because of* the insufficiency of angular 
momentum. It follows that only a part of the uplifted matter from the star 
can finally emerge to become the ring and the rest after having transferred 
a part at least of their angular momentum to the particles in the ring must 
either escape into infinity or fall back to the star. This idea advanced in 
1973 was developed recently into a quantitative theory (Huang 1976) by the use 
of laws of conservation of energy and angular momentum. The theory shows that 
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if a rotating star ejects matter, ring formation is a consequence consistent 
with the lavs of conservation of energy ^nd angular momentum. 

Once formed, the ring is not an isolated system apart from the gravita- 
tional attraction of the central star, because the latter is continually, if 
not continuously, ejecting new mass, which could Interfere with the gas al- 
ready in the ring. Thus, the ring structure will be modified or even destroyed 
by its interaction with the newly ejected matter, depending upon the degree 
of violence of ejection. The present paper is to present a quantitative de- 
scription of the modification of the ring structure by the newly ejected matter. 
It is then followed by a demonstration of such modifications in two stars-- 
P Mon and 7^ Aqr — which have so far defied satisfactory interpretation 
based on a ring with constant semi-major axis and constant eccentricity. 

II. MODIFICATION OF THE RING STRUCTURE 

In what follows we will use dimensionless variables by adopting the break- 
up rotational velocity at the equator of the star as the unit of velocity and 
its equatorial radius as the unit of length. Hence the angular momentum per 
unit mass of the matter moving at the break-up velocity in the equator will 
also be the unit of angular momentum per unit mass. In this system of units 
we denote by and respectively the tangential and the radial component 
of ejected velocity of the matter at the equator, and by m, k and e respectively 
the total mass, the semi-major axis and the eccentricity of the gaseous ring 
around the star. All of the quantities m, k,. and e are considered to change 
with time either abruptly or slowly as the case may be. Finally it may be 
noted that the ring must have a lateral extent. So k should be considered as 
an average value. 

Before we formulate the relations that govern the modification of the 
ring structure by the newly ejected matter from the star we must first state 
that not all ejected matter will have an elfect on the ring, because a part 
of this matter falls back to the star before it reaches the ring. After the 



ejected matter has mixed with the gas in the ring, the mass will be further 

dissipated as some particles will fall back to the star and others escape into 

infinity. Thus if we let be the ejected mass, A m the increase of the 

mass in the ring and A m. the mass dissipated (i.e. either falling back to 

d 

the star or escaping into infinity), we hav Am," A m+ Am.. It may be 

J a 

pointed out that m can be positive, zero, or negative, because there is 
no a priori reason that A m^ must be greater than, equal to, or less than A 
We now assume that the ejection of A ® j star affects the ring 

structure such that its total mass is changed from m to m+ A m, its semi-major 

, • 

axis from k to k+ A h and its eccentricity from e to e+ A e. Since A mj is 

the cause of the modification we may use m^ as the independent variable. If 

furthermore the ejection is continuous and at a constant rate, measures also 

the time if the rate of ejection Is known. 

In our unit system the specific angular momentum of A m^ is x^. Let us 

assume the specific angular momentum of the dissipated mass, namely Am^, to 

be h .. A consideration of angular momentum before and after the star has 
d 

ejected A gives 






I 


( 1 ) 


or 




( 2 ) 


when it is expressed in a differential form, where we have set 
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which represents the angular momentum input to the ring per unit ejected masst 
Similarly a consideration of the dyaiamical energy before and after the 
ejection of ^ by the star gives 




(A) 


where all terms except the last one in the right hand side are self-evident. 
This last term includes the dynamical energies (1) dissipated into heat, (2) 
carried away by the falling and escaping mass ^ m^, and (3) obtained by the 
work done by the heat energy of the ejected matter in the process. Under the 
usual condition b is a positive number, because in the process of mixing of 
the ejected streams with the gaseous ring the dissipation and the escape of 
dynamical energy must be more effective than the conversion of heat into dynam 
ical energy. We may again write equation (4) in a differential form as 


where 


0 

c, = ■ 


(5) 


( 6 ) 


which represents the dynamical energy that the ring dissipates for each unit 

of ejected mass. Thus we see that the modification of the ring structure by 

the newly ejected matter is controlled by two parameters, the angular momentum 

input factor and the energy dissipation factor . 

If we denote by m , k and e the Initial values of m, k and e respectively 
o o o 

of the ring before the ejection has started, i.e. when m^* 0, integration of 
equations (2) and (5) yields 
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■m * ™o £e ’"j 


and 


7n % . ^ ^ 


whe re 


E. 


(7) 


( 8 ) 


(9) 


Here A and C are respectively values of Aj^ and Cj^ averaged over the entire 
process in which a total mass of has been ejected from the star. 

In general we may write 


w = /rnj 


( 10 ) 


where ^ is the fraction of ejected mass that enters into the ring, Hepce 
^ is at most equal to 1 but can be negative if the ejected matter induces 
the ring to dissipate more mass than is acquired. If we further write 



( 11 ) 


k and e of equations (7) and (8) may be expressed in terms of ^ which replaces 
m^ as the independent variable in the following discussion. 

In order to show the effect of ejected mass on the ring structure most 

(i^) (lij 

comprehensively we eliminate I; and m, from equations (7), (8)^^and to give 
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C as a function of p , vlth e as a parameter: 




( 12 '> 


Thus C(^; e) constitutes a family of curves on the (pt* C) plane, each curve 
being labelled by a value of e. In figure 1 we illustrate this family of 
curves for the case 0 (i.e. the mass of the gaseous ring remaining con- 
stant) in the main figure and for the case ^ 1 (i.e, complete merging of 

the ejected mass with the ring) in the inset. For the case ^ " 0, we have 

dm,/dm<= 1, Since we do expect that the dissipated mass carries away some 

d J 

angular momentum, we have Illustrated this case in the figure with A ■ 0.6. 

In the case 3=1, dm./dm-* 0, we have illustrated this case in the inset with 

d j 

A = 1. Actually, as will be seen later, the value of A is immaterial to the 
general conclusion that can be reached from the figure. In both illustrations 
we have adopted k^= 3 and e^= 0.25. Observatlonally the case of constant 
mass ( ^ *0) is closer to the reality than the case of complete merging 
-C / =1)> at least during the time in which the total emission strength of 
a line remains constant (e.g. Struve 1931, McLaughlin 1961). 

The single-parameter family of curves C(^; e) as labelled by parameter 
e is divided by the critical curve C(^ ; e^) into two groups. Above the critical 
curve are curves in the family with e^e and below it are curves in the 
family with e^e^. In the upper region all curves of e approaches + OC as 
^ 0 while in the lower region all curves of e^e^ approaches - OC as 

^ 0. The significance of this family of curves can be seen when we pick 

a value of C and trace the horizontal line (i.e. C = constant) in the figure. 
Since the ring starts with e = e^ , at p “ 0, the value of e decreases or in- 
c reases^ as 


below 


p in'creases at any point, depending upon whether C is above or 
C(p, bet us take the case of C “ 0. Following C • 0 horizontally 
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in figure 1, we intercept first curve C(^; 0,2) at ^ ■ 0,033, then curve 
C<p; 0.1) at ^ " 0.077 9 nd finally curave C(^ ; 0) at ^ ■ 0.09. This shows 
that as the ejection continues, the eccentricity decreases, until it becomes 
circular. Once the ring is circular, the level of dissipation as given by C 
cannot be maintained, because with e “ 0 the kinetic energy of the gas in the 
ring is already completely associated with the angular momentum and cannot be 
dissipated unless the angular momentum input A assumes a different value. 
Therefore the region above C(^, 0) is forbidden. Riysically before the ring 
becomes circular, the energy dissipation and the angular momentum input will 
adjust themselves, so it is not realistic to envisage the course of ring 
modiilcation under the condition of constant C and A at that time. One can 
easily see that when the eccentricity is close to zero, either the dissipation 
has to dei rease or the angular momentum input has to decrease or both. This 
means that C and/or A will decrease as the eccentricity approaches zero. 

On the other hand if the C is small, say C » -0.2, we find in the similar 
way as just described that the eccentricity of the ring increases as in- 
creases along the line C “ -0.2, But if C lies only a little above the point 
given by C(0, e^) , we find that the eccentricity of the ring decreases first, 
then increases, crosses the C(^, e^) curve and finally Increases as in- 
creases, because C(^, e^) is not exactly horizontal but is slightly slanted 
upward with increasing ^ . But the change in e with ^ in this instance is 

very small at any value of p . 

From these we can see clearly that for the case 3*0 and A " 0.6 that 
is shown in the figure any C^-0.1 tends to decrease the ring eccentricity and 
any C <. -0.12 tends to increase the eccentricity. In general large dissipation 
decreases e and small dissipation increases e. 

Let us consider now the effect of the angular momentum input factor A. 

For samll values of ^ we derive from equation (12) that: 
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c-^ 




(13) 


Thus for each curve in figure 1 (given e^, £^, e, and }() a change of A and/or 
^ simply translates It vertically, If we confine ourselves to small ^ . 

For a constant ^ , an Increase In A corresponds to a translation of all 
curves downward and a decrease in A a translation upward. Otherwise the rela* 
tlve positions of different curves are not changed. This means that an in- 
crease in the angular momentum Input has the same effect on the eccentricity 
of the ring as an increase in the dissipation. Similarly for a given value of 
A, a change in ^ simply translates the entire family of curves upward or 
downward depending whether ^ increases or decreases if ^ remains small. 

This means that an increase in ^ has the same effect on the eccentricity of the 
ring as a decrease in the dissipation. 

So far we have limited our discussion only to cases of small ^ , which 

corresponds to mild ejections that slowly modify the ring structure. Wliat 
would happen if the ejection is so violent that the ejected mass is much greater 
than the total mass in the ring? Figure 1 provides us with an answer to this 
question. If the energy dissipation is high, the eccentricity will be abruptly 
reduced to zero. But if the energy dissipation is small, a violent ejection 
abruptly increases e beyond 1 and the ring is completely destroyed. 

The variation of k with 3 studied in the s^ way as that of e. 

Thus we may express from equations (8), (10) arid (11) 


C 




f 


i±<l. 1 

I L 


(14) 


whiclt forms another family of curves on the (jJ-C) plane with k as the labelling 
parameter. It can be easily seen that this family of curves is divided Into 
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two groups by a horizontal line 



( 15 ) 


which Is a critical member of the family. All C(p; k) curves corresponding 

to k>k He below this line and all C(&; k) curves corresponding to k<k 
o r o 

lie above that line. All these curves In the family. whether k>k or k^ k 

/ o ^ o» 

are hyperbolas with the two lines 




0 




.(16) 


as asymptotes. Hence for k<k^, the hyperbolas stack upward as k decreases 
from k and for k^k , the ' 'perbolas stack downward as k increases from k . 
Thus as ^ increases^k increases indefinitely when C 4, ^ / (2k^) , and decreases 
indefinitely when C / ?(/(2k^). In the special case 0, k Increases or de- 

creases according to whether C is negative or positive. Thus we see that the 
two families of curve C(^ ; e) and C(^ ; k) are very similar. Each family Is 
divided into two groups by a critical curve In its own family, namely C(^; e^) 
and C(^ ; and k^ being the initial values of e and k of the ring. 

By copibining the family of curves C(^; e) , such as In figure 1, with that 
of curves C(^; k), one can obtain at a glance the changes in both k and e at 
any j; for different values of C without actually computing them. In order 
to do so we may simply draw the two critical curves C(^, e^) and C(^; k^) in 
the ( -C) diagram. Then the region In the (^-C) diagram above the upper 
one of these two critical curves corresponds to C and ^ values that make 
both k and e decrease as ^ increases. The region below the lower one of the 
two critical curves corresponds to C and ^ values that make both k and e in- 
crease as B Increases. The region between the two curves Includes those C 
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and p values that will make either k increase and e decrease or k decrease 
and e increase, depending upon whether 0(0 ; k ) is above or below C(B , e ). 

Y o r ° 

In the case of A ■ 0.6, 3 " figure 1 shows clearly C(^ , k^) which corres- 

ponds to the horizontal line, C ■* 0, is above C(^ ; e^). So we have in the 
intermediate region increasing k and decreasing e as ^ increaseb; while in 
the case of A ■ 1, ^ ■ 1, the inset in figure 1 shows that C(^ ; k^) which 
corresponds to the horizontal line C below C(p, e^). So the 

intermediate region on the ( ^ -C) plane includes C and ^ values that make 
k decrease and e increase as the star ejects mass. 

In this way we have seen how the modification of the ring structure is 
controlled by three parameters, C, A and ^ and found that all combinations of 
the changes of k and e are possible. Hence we cannot rule out any of the 
combinations purely from the present consideration. > 

ITT. OBSERVATIONAL EVIDENCE OF ORBITAL MODIFICATION 
Before we interpret the observational results i.' the light of ring mod- 
ification just presented, we may first summarize what has been done with the 


rotating eccentric ring without orbital modification. First we have devel- 
oped a quantitative theory of the eccentric ring (of constant e and k) revolving 

around Be stars with a slow drift of Ics line of apsides (Huang 1973a). Since 
from 

it takes^several years to more than a decade for the line of apsides to complete 
a revolution, the motion of emitting atoms in the r^ng may be regarded as 
following Kepler's laws of motion during the 'exposure time that each spectro- 
gram was taken. Thus from tl V/R variation of the emission line and radial 
velocity curves of the emission edgec and the central absorption of the line. 


we are able| to determine the eccentricity and the semi-major axis of the ring. 

a 

Altogether four stars have been studied in two previous papers (Huang 1973,^«ml' 
Albert and Huang 1974) from this point of view. From the values K ^(In km s ^), 
which is one half of the separation between red and violet edges of the emission 
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line, and e, the eccentricity of the ring, determined from obaervatlona In 

those two papers we have recalculated k values (which Is denoted by a^/R In 

.1 

the previous papers) by adopting uniformly ''10 km a (Slettebak 1966) as the 

-1 

break-up rotational velocity, V, for B3 stars and S62 km s as that for 
B2 stars and 1.8 as the mass-radius ratio (In the solar unit) for B3 stars 
and 2.0 as that for B2 stars. We have estimated the value of sin 1, where 
i denotes the inclination of the plane of the ring which Is assumed to be 
Identical to the plane of the equator, from the observed rotational velo- 
city V sin 1 and the adopted' break-up velocity V. Since according to our 
Interpretation, rings can be formed even when the equatorial rotational 
velocity is less than V, we may have underestimated sin 1, and therefore 
k. On the other hand according to our theory 2K^ should be the difference 
between the outermost edges of emission lines. But in the actual measure- 
ment, the edge is determined by the steepest slope of the density change on 
the spectrogram, which is way inside the actual edge. Consequently is 

greatly and systematically underestimated. An undereat Imate for K causes 

e 

an over-estimate of k. 

Table 1 gives the result of this uniform calculation where P (in 

a 

years) is the period of apsidal motion. Only tl»ree stars are listed here 
because the fourth, p Mon, is re- ir » erpreted later in this paper. From 
the figures in the previous two papers it can be seen that these three 
stars show good agreement between observations and calculated results for 
both the V/R variations of the emission line and the velocity curves of the 
emission edges and central absorption. However in ^ ^k)n, we have found 
that the calculated curves do not agree well with the observed result ob- 
tained by Cowley and Gugula (1973). The data of another star, “JT Aquaril, 
obtained by KcLauglilin (1962), completely defy the Interpretation based on 
the rotating eccentric model of constant k and e. So no attempt was even made 
by us to underst and t.iis star simply on the basis of this moilel . Indeed the 
difficulty found in these two stars prompted us to think 
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chat some modification of Che rotating eccentric iiKidel es has been presented 
in the previous section is necessary. In light of the modification of the ring 
we can now proceed to exuiine the observed results for ^ Mon and XAqr. 

IV, IMTERPRETATIOK OF A* HONOCEROIIS 
The major difficulty we have found in Che previous invest:' ition of this 
star is chat no ring with fixed values of k and e can satisfy Che observation. 
But if we are prepared to consider the ring as having undergone a few changes 
in k and e during the period from 1930 .o 1970, the eccentric rotating ri^g 
can explain the observational data much better than that with constant k and 
e. The new interpretation is shown in figure 2 which may be compared with 
figure 3 in the previous paper (Albert and Huang 1974). The improvement is 
obvious. 

! gure 2 where the observed points were taken directly from Cowley 
and Gugula (1974) we have interpreted the gaseous ring as has gone through 
three noticeable changes in the time of observation from 1930 to 1970. Three 
vertical lines show the approximate times of these changes. From the observa- 
tional data we have found that each change occurs in a relatively short inter- 

val so that the calculation (Huang 1973a) based on the eccentric rotating ring 
with constant k and e is applicable to the period between any two consecutive 

changes, but we allow for some differences in the values of k and e in differ- 

ent periods. This is equivalent say that the ejection that is strong enough 
to modify the ring structure is sporadic rather than continuous. Or more 
specifically the time interval that the ring structure is being modified by 
such an ejection is short compared with the period of apsidal motion. How- 
ever it may or may not be short compared with th? orbital period. Previously 
we have suggested that disappearance of emission lines at certain times in 
some Be stars is also due to sporadic ejections which destroy the ring com- 
pletely (Huang 1973b). It goes without saying tiiat ejections iliat modify tlie 
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ring must be mild compared with those that destroy the ring. 

Ir order to present the results o£»our analysis, we denote by A, B, 

C, and D the four periods separated by three vertical lines that correspond 
to the respective times of ring modification. The values of P^, k and e 
determined in each epoch are based mainly on the V/R variation. Solid curvea - 

representing the V/R variation were first plotted from the determined value^ 

of e. The velocity variations are not as clear-cut as the V/R variation, as 
we can easily see from figure 2. This is understandable if we consider the 
fact that the value of V/R is determined by the ratio of the number of onit- 

ting atoms in the entire approaching portion of the ring to that in the en- 
tire receding portion, while the three radial velocities are determined by, 

only small portions of the matter in the ring. Thus any turbulence can 

greatly distort the radial velocity cf particles in the small portion of the 
ring, but its effect on the V/R variation is somewhat reduced because large 
portions of the ring are involved. This sinply reflects the statistical 
truth that a small (large) sample corresponds to a large (small) fluctuation. 

It appears that observed radial velocities of the central absorption over 
the entire period of observation can be best fitted to a straight line slightly 
tilted downward. It means that the radial velocity has decreased a few km s ^ 
in the forty more years. Since this represents, according to the theory^ 

the radial velocity of the star, such a decrease might represent 
the orbital motion of the star around the distant companion, if the decrease 
should be intrinsic. At present we are not sure whether this is so. 

Once radial velocity corresponding to the central absorption is deter- 
mined, the K value in each epoch is obtained by finding two horizontal lines 
e 

placed symmetrically above and below the central-absorption one so that the 
two fit the observed data respectively of the red and violet emission edges 
best. The predicted velocity curves for the red and emission edges and central 
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absorption can then be calculated from and e and are shown as the solid 

cup?es in the figure. • 

* 

In table 2 the values, P , k and e thus determined are given for each 

of these four periods. The eccentricity of the ring in epoch C te^not w l A 

unceiHa'v'. 

determined. So 4^ i^ |lven inside parentheses^^ It appears that from epoch 
A to epoch B, both the eccentricity and the semi-major axis have increased. 
But from epoch B to epoch C, the eccentricity ' is e o d w e ad and ^finally be- 
comes 0 at epoch D, while the semi-major axis remains practically cpnstant 
during all these three epochs. 


V. T AQUARII , 

The complexity of the V/R variation and the accompanying radial vel- 
ocity curves of emission edges and the central absorption of 7T Aquarii 
has induced to propose the ring m'^ diflcation. Indeed once the idea of 
ring modifir 'i -ion is accepted, we have found an easy explanation of the seem- 
ingly irregular data obtained by McLaughlin (1962) in the period of half a 
century. Our interpretatioi* of these data is illustrated in figure 3 where 

circles represent the observed . ormal points given by McLaughlin in Table 2 

Q ^ 

of his paper. Following what has been done in the previous section for u 
Mon, we have divided the fifty years into thirteen separate epochs denoted 
by A, I, C, . . . M. This shows that “JT Aqr is an active star and ejections 
of matter that are strong enough to modify the ring structure occurred quite 
often. 

In order to determine k from K , we must know sin i (where i is the in- 

e 

clination of the plane of the ring) and the mass-radius ratio in the solar 
unit. According to McLaughlin, 7T Aqr is a Bl star. So we take the mass 
radius ratio to be 2.1. In the Appendix we have estimated the observed 
rotational velocity V sin 1 to be 480 km s ^ . Combined with break-up rota- 
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tlonal velocity V equal to 590 km s (Slettebak 1966) It gives a rough es' 

timate of sin i ■ .814. Table 3 gives th£ values P , e, K and k detemined 

a e 

from the data for different epochs. Observed points of most epochs and the 
corresponding c,urves and lines are shown in figure 3. A few >t{pochs (A» H* 

I, J, K) are not illustrated because there is either no emission or 

no V/R variations (i.e. t zero eccentricity). 

VI. DISCUSSION 

With the understanding that the ring structure may be modified from 
time to time, we can now see that stars whose V/r variations have been 
studied intensively all show this behavior. The only difference is the 
frequency of these modifications. For stars, like HD 20336 and 25 Ori, the 
regular periodic variation of V/R and of the radial velocities of emission 
edges and central absorption can last a few cycles without showing any 
sign of significant modifications. But Aqr which shows frequent changes 
in its emission behavior, is, according to our interpretation, more active 
in its ejection of matter than other stars studied. 

The fact that the data derived from Be stars can all be understood so 
far in terms of a ring of constant k and e in each Interval between two con- 
secutive modifications shows that the disturbances that modify the ring struc 
ture last only short times compared with the quiescent periods even for J:he 
most active stars observed. That is why we may vii(^ the ejections that take 
place on the surface of Be stars as sporadic.* From observations it appears 
that the time interval between two sporadic ejections is on the average of 
the order of a few years to a few decades. Whether there are Be stars in 
which sporadic ejections take place so often as to become practically con- 
tinuous, we have no idea. If such stars should exist, the V/r variations 
would become erratic and no analysis of the kind used here can be applied. 


The observed data also led us to believe that one and the same star 



could be quiescent for a long time and then become active In its ejection 
of matter later. At present it is not* clear whether there is any regular 
cycle of quiescence and activity for each star. 

While the apsldal motion explains fairly well the V/R variation and 
velocity curves of emission edges and central absorption, we doubt, after the 
present study, that the apsldal motion is due solely to the equatorial bulge 
of the star, as it has been thought to be so far (Struve 1931, Johnson 1958). 
The reason for our doubt arises from the argument that if the apsldal motion 
should b^due to the equatorial bulge, its rate of advance would depend 
critically on the semi-major. axis of the ring (e.g. Sterne 1939). As a 
result the larger the orbit, the less would be the rate of apsldal motion 
and the longer would be the period P . Such a trend has been found neither 
in p Mon nor in *71 Aqr If we examine tables 2 and 3 carefully. If this 
trend were not obscured by errors In the determination, we would suggest 
that the apsidal advance may be influenced by other factors in addition to 
the gravitational effect of the equatorial bulge of the star. 

The present investigation has been supported by a grant from National 
Aeronautics and Space Administration. 


APPENDIX 


Rotational Velocity of TIT Aquaril 

According to McLaughlin (1962), 7T Aqr is a Bl star In whose spectrum 
Mg II X 4481 "can be measured on plates of the best contrast." But Rlng- 
uelet and Machado (1974) have classified the stars as BO and consequently 
considered He I ^ 4471 as an asymmetrical line instead of a blend of this 
line with Mg II ^ 4481. They have reproduced the tracing of this blend 
complete with the wavelength scale in their paper in order to show the asym- 
metry of the line. However when we have followed McLaughlin's interpretation, 
we can clearly draw from their tracing two rotationally broadened symmetric 
lines each corresponding roughly to an observed rotational velocity (V sin i) 
of about 480 km s ^ and with a separation that corresponds to that between 
He I ^ 4471 and Mg II 4481 . Indeed if one should consider the entire 
blend as one single line, the rotational velocity would amount to more 
than 700 km s ^ which obviously is not reasonable. 


o IT Kjf i'< /iXj ife -• 
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Legends 

Fig. 1. — Energy dissipation factor C(^; e) as a function of ^ with e as 

the labelling parameter. The critical member of this family of curves Is 

C(A; e ) where e is the initial value of e when the newly ejected matter 
VO o 

starts to modify the ring structure. 

Fig. 2. — Interpretation of the observed results (1930-72) of ^^Mon on the 
basis of ring modifications. The points are observed values taken from 
Cowley and Gugula’s (1973) paper while the curves and lines represent our 
interpretation given in the text. 

Fig. 3. --Interpretation of the observed results (1911-61) of “JC Aqr on the 
basis of ring modifications. The open circles are observed values taken from 
McLaughlin's (1962) paper, while the curves and lines are our interpretation 
given in the text. 



Table 1 


Characteristics of Gaseous Rings of Three Stars 


Star 

epoch 

^(year) 

K (km s 
e 

e 

k 

105 Tfluri 

1930-52 

11.5 

132 

.2 

3.7 

HD 20336 

1916-32 

4.5 

198 

.25 

4.1 

25 Orionis 

1914-25 

4.8 

168 

.26 

4.7 
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Table 2 

Evolution of the Gaseous .Ring Around Mon 


Designation 

Epoch 

g^(year) 

K (km s 
e 

e 

k 

A 

1930-54 

12.8 

200 

• 3 

4.9 

B 

54-64 

11.2 

182 

-4 

6.4 

C 

64-68 

n.d. 

171 

n.d. 

(6.4) 

0 

68-72 


165 

0 

6.5 


n.d.: not determinable 


Table 3 


Evolution of the Gaseous Ring Around t Aqr 


Designation 

epoch 

^(years) 

Kg (tan s“S 

e 

k 

A 

1911-23 

— 

205 

0 

6.8 

B 

24-28 

5.1 

218 

.22 

6,4 

C 

29-32 

3.3 

225 

.4 

6.8 

D 

33-36 

3.6 

241 

.35 

5,6 

E 

37-38 

no emission 




F 

39-41 

— 

261 

0 

4.2 

G 

41-44 

3.3 

214 

.25 

6.7 

• H 

44-45 

no emission 




I* 

47-49 

— 

351 

0 

2.3 

J 

50 

no emission 




K 

51-54 

— 

250 

0 

4.6 

L 

55-57 

4.0 

212 

.2 

6.7 

M 

58-61 

5.5 

212 

.25 

6.8 


’^his epoch may have witnessed two changes In k. But^ not enough observational 
points to justify a more refined determination. The eccentricity has remained 
to be zero during the entire epoch. 
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kin/MC. 




